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Abstract — In  this  paper  we  study  the  problem  of  coverage  in 
heterogeneous  planar  sensor  networks.  Coverage  as  a  perfor¬ 
mance  metric,  quantifies  the  quality  of  monitoring  provided  hy 
the  sensor  network.  We  formulate  the  problem  of  coverage  as  a 
set  intersection  problem  arising  in  Integral  Geometry,  and  derive 
analytical  expressions  for  stochastic  coverage.  Our  formulation 
allows  us  to  consider  a  heterogeneous  sensing  model,  where 
sensors  need  not  have  an  identical  sensing  capability.  In  addition, 
our  approach  is  applicable  to  scenarios  where  the  sensing  area 
of  each  sensor  has  arbitrary  shape  and  sensors  are  deployed 
according  to  any  distribution.  We  present  analytical  expressions 
only  for  convex  sensing  areas,  however,  our  results  can  be 
generalized  to  non-convex  areas.  The  validity  of  our  expressions 
is  verified  by  extensive  simulations. 

I.  Introduction 

Sensor  networks  are  projected  to  have  a  significant  impact 
into  our  everyday  lives,  with  applications  to  environmental 
monitoring,  home  health  care,  disaster  relief  operations,  and 
ambient  monitoring  [1].  One  of  the  primary  tasks  of  sensor 
networks  is  the  collective  monitoring  of  a  field  of  interest. 
Sensors  may  monitor  physical  properties  such  as  temperature, 
humidity,  air  quality,  or  track  the  motion  of  objects  moving 
within  the  field  of  interest.  In  order  for  the  sensor  network 
to  sufficiently  monitor  the  entire  field  of  interest,  one  needs 
to  ensure  that  every  point  of  the  field  is  covered  by  at  least 
one  sensor.  Eurthermore,  to  provide  the  desired  accuracy  and 
robustness  against  node  failures,  many  applications  require 
that  each  point  of  the  field  of  interest  is  sensed  by  more 
than  one  sensor.  Hence,  the  problem  of  node  deployment  for 
the  purpose  of  sensing  can  be  viewed  as  a  coverage  problem, 
defined  below. 

The  coverage  problem  is  to  quantify  how  well  is  the  field  of 
interest  sensed  by  the  deployment  of  the  sensor  network.  The 
coverage  problem  can  be  studied  under  different  objectives  and 
constraints  imposed  by  the  applications  such  as,  worst-case 
coverage  [10],  deterministic  coverage  [10],  [12]  or  stochastic 
coverage  [8],  [10],  [12],  [15],  [20].  The  worst-case  coverage 
problem  quantifies  coverage  based  on  the  parts  of  the  field  of 
interest  that  exhibit  the  lowest  observability  from  the  sensors 
[10],  and  is  relevant  in  applications  where  a  desired  threshold 
number  of  sensors  need  to  observe  the  field  of  interest. 
The  deterministic  coverage  problem  [10],  [12]  quantifies  the 
coverage  achieved  by  deploying  sensors  in  a  deterministic 
way,  and  is  relevant  in  applications  where  one  can  select 
the  positions  where  the  sensors  are  placed.  The  stochastic 
coverage  problem  [8],  [10],  [12],  [15],  [20],  on  the  other  hand, 
quantifies  the  coverage  achieved  when  sensors  are  deployed 


according  to  a  distribution,  and  is  relevant  in  applications 
where  the  sensors’  positions  cannot  be  selected  a  priori. 

In  this  paper,  we  analyze  the  following  stochastic  coverage 
problem.  Given  a  planar  field  of  interest  and  N  sensors 
deployed  according  to  a  known  distribution,  compute  the 
fraction  of  the  field  of  interest  that  is  covered  by  at  least 
k  sensors  (k  >  1).  The  problem  can  also  be  rephrased  as, 
given  a  field  of  interest  and  a  sensor  distribution,  how  many 
sensors  must  be  deployed  in  order  for  every  point  in  the  field 
of  interest  to  be  covered  by  at  least  k  sensors  with  a  probability 
p  (k-coverage  problem)  [20]. 

In  this  paper  we  make  the  following  contributions.  We 
formulate  the  problem  of  coverage  in  sensor  networks  as 
a  set  intersection  problem.  We  use  results  from  integral 
geometry  to  derive  analytical  expressions  quantifying  the 
coverage  achieved  by  stochastic  deployment  of  sensors  into 
a  planar  field  of  interest.  Compared  to  previous  analytical 
results  [8],  [12],  [20],  our  formulation  allows  us  to  consider 
a  heterogeneous  sensing  model,  where  sensors  need  not  have 
an  identical  sensing  capability.  In  addition,  our  approach  is 
applicable  to  scenarios  where  the  sensing  area  of  a  sensor 
is  not  an  ideal  circle,  but  has  any  arbitrary  shape.  To  the 
best  of  our  knowledge,  only  [15]  considers  a  heterogeneous 
sensing  model,  though  only  incorporating  the  mean  value  of 
the  sensing  range  in  the  coverage  computation.  In  addition,  the 
formulation  in  [15]  considers  only  uniformly  deployed  sensors. 
In  our  approach,  sensors  can  be  deployed  according  to  any 
distribution.  We  provide  formulas  for  k-coverage  in  the  case 
of  heterogeneous  sensing  areas,  as  well  as  the  simplified  forms 
in  the  case  of  identical  sensing  areas,  and  give  an  example  for 
the  computation  of  the  number  of  sensors  required  to  cover 
a  field  of  interest  with  a  pre-specified  probability.  Einally,  we 
validate  our  theoretical  expressions  via  simulations  and  show 
an  exact  match  between  simulation  and  theory. 

The  rest  of  the  paper  is  organized  as  follows.  In  Section 
III  we  formulate  the  coverage  problem  as  a  set  intersection 
problem.  In  Section  IV  we  derive  analytical  expressions  for 
coverage.  In  Section  V,  we  validate  our  theoretical  results  via 
simulation.  Section  VI  presents  our  conclusions. 

II.  Related  Work 

In  this  section  we  describe  related  work  to  the  coverage 
problem  in  wireless  sensor  networks.  The  coverage  problem 
can  be  classified  under  different  objectives  and  metrics.  The 
different  approaches  to  the  coverage  problem  are,  deterministic 
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or  stochastic  sensor  deployment,  homogeneous  or  hetero¬ 
geneous  sensing  area,  additional  design  constraints  such  as 
energy  efficiency,  minimum  number  of  sensors  that  need  to 
be  deployed,  or  network  connectivity.  Based  on  the  objective, 
the  coverage  problem  formulation  varies  to  reflect  the  different 
assumptions  and  objectives. 

In  [9],  the  authors  study  the  problem  of  deterministic  node 
placement  in  order  to  achieve  connected  coverage,  that  is, 
sense  the  field  of  interest  with  the  minimum  number  of  sen¬ 
sors,  while  keeping  the  sensor  network  connected.  The  authors 
assume  that  the  sensing  area  of  each  sensor  follows  the  unit 
disk  model  and  consider  sensors  with  identical  sensing  areas. 
The  problem  of  connected  coverage  has  also  been  recently 
studied  in  [20].  The  authors  provide  a  geometric  analysis 
that  relates  coverage  to  connectivity  and  define  the  necessary 
conditions  for  a  network  covering  a  field  of  interest  to  be 
connected.  The  conditions  for  coverage  and  connectivity  are 
derived  based  on  the  assumptions  that  the  sensing  area  of  each 
node  is  identical  and  circular,  and  the  location  of  the  nodes 
is  known.  The  authors  extend  their  algorithms  for  the  case  of 
probabilistic  deployment,  and  also  relax  their  assumptions  to 
non-unit  disk  sensing  areas,  by  approximating  the  real  sensing 
area  with  the  biggest  possible  circular  area  included  in  the  real 
sensing  area. 

In  [16]  the  authors  study  the  problem  of  deterministic 
coverage  under  the  additional  constraint  that  each  sensor  must 
have  at  least  k  neighbors.  They  propose  a  deployment  strategy 
that  would  maximize  the  coverage  while  the  degree  of  each 
node  is  guaranteed  to  be  at  least  k,  under  the  assumption  that 
the  sensing  range  of  the  sensors  is  isotropic. 

In  [13],  the  authors  study  the  problem  of  coverage,  as  a 
path  exposure  problem.  Using  a  generic  sensing  model  and 
an  arbitrary  sensor  distribution,  they  propose  a  systematic 
method  for  discovering  the  minimum  exposure  path,  that  is  the 
path  along  which  the  network  exhibits  the  minimum  integral 
observability^.  Authors  in  [10],  investigate  the  problem  of 
best-  and  worst-case  coverage.  In  their  formulation  of  the 
coverage  problem,  given  the  location  of  the  sensors  and  a 
generic  sensing  model  where  the  sensing  ability  of  each 
sensor  diminishes  with  distance,  the  authors  use  Voronoi 
diagrams  and  Delaunay  triangulation  to  compute  the  path  that 
maximizes  the  smallest  observability  (best  coverage)  and  the 
path  that  minimizes  the  observability  by  all  sensors  (worst 
coverage).  In  [11],  the  authors  provide  a  decentralized  and 
localized  algorithm  for  calculating  the  best  coverage. 

Authors  in  [12],  study  the  problem  of  stochastic  coverage 
in  large  scale  sensor  networks.  For  a  randomly  distributed 
sensor  network,  the  authors  provide  the  fraction  of  the  field 
of  interest  covered  by  k  sensors,  the  fraction  of  nodes  that 
can  be  removed  without  reducing  the  covered  area  as  well 
as  the  ability  of  the  network  to  detect  moving  objects.  The 
results  presented  in  [12]  hold  only  for  randomly  (uniformly) 
deployed  networks  and  under  the  assumption  that  the  sensing 

*The  integral  observability  is  defined  as  the  aggregate  of  the  time  that  a 
target  was  observable  by  sensors  while  traversing  a  sensor  network. 


area  of  each  sensor  is  identical.  Furthermore,  the  analysis  in 
[12]  suffers  from  the  border  effects  problem,  illustrated  in  [2], 
[3].  The  results  hold  asymptotically  under  the  assumption  that 
the  field  of  interest  expands  infinitely  in  the  plane,  while  the 
density  of  the  sensor  deployment  remains  constant. 

In  [15],  the  authors  study  the  stochastic  coverage  problem 
in  ad  hoc  networks  in  the  presence  of  channel  randomness.  For 
a  randomly  deployed  sensor  network,  the  authors  analyze  the 
effects  of  shadowing  and  fading  to  the  connectivity  and  cov¬ 
erage.  They  show  that  the  in  the  case  of  channel  randomness, 
the  coverage  problem  can  still  be  modeled  after  the  Spatial 
Poisson  distribution,  by  using  expected  size  of  the  sensing 
area  of  sensors.  While  the  results  in  [15]  are  applicable  to 
heterogeneous  sensor  networks,  they  hold  only  for  randomly 
deployed  networks,  and  are  impacted  from  the  border  effects 
problem  [2],  [3],  as  noted  by  the  authors  [15]. 

Compared  to  previous  work  that  derives  analytical  coverage 
expressions  [12],  [15],  [16],  our  formulation  allows  us  to 
consider  a  network  model  where,  (a)  sensors  can  be  deployed 
according  to  any  distribution,  (b)  sensors  can  have  a  sensing 
area  of  any  arbitrary  shape,  (c)  sensors  can  have  heterogeneous 
sensing  areas. 

III.  Problem  formulation  &  Background 

In  this  section,  we  formulate  the  problem  of  coverage  in 
heterogeneous  sensor  networks  as  a  set  intersection  problem 
arising  in  Integral  Geometry  [7],  [14],  [17]-[19]  and  provide 
relevant  background  for  the  set  intersection  problem. 

A.  Problem  Formulation 

We  formulate  the  problem  of  stochastic  coverage  as  follows. 
Let  ^0  denote  the  planar  field  of  interest  we  want  to  monitor, 
with  area  Fq  and  perimeter  Lq.  Assume  that  N  sensors  with 
sensor  Si  having  a  sensing  area  Ai,  {i  =  l...A^),  are 
deployed  according  to  a  distribution  K{Ao)  and  in  such  a  way 
that  they  sense  some  part  of  the  field  of  interest^.  Let  Fi,  Li 
denote  the  size  and  the  perimeter  of  the  sensing  area  Ai  of 
each  sensor  Si,  respectively.  We  want  to  calculate  the  fraction 
of  Aq  that  is  sensed  by  at  least  k  sensors,  i.e.  the  fraction  that 
is  k-covered  (k  >  1).  This  problem  is  equivalent  to  computing 
the  probability  that  a  randomly  selected  point  P  G  is 
sensed  by  at  least  k  sensors.  We  map  this  coverage  problem 
to  the  following  set  intersection  problem.  In  our  formulation, 
a  set  S  is  defined  as  a  collection  of  points  in  the  plane,  and 
for  the  coverage  problem  the  sets  are  closed  regions. 

Let  Sq  be  a  fixed  closed  set  defined  as  a  collection  of 
points  in  the  plane,  and  let  Fq  and  Lq  denote  the  area  and 
perimeter  of  Sq.  Let  N  closed  sets  Si  (i  =  1 ...  N)  of  size  Fi 
and  perimeter  Li  be  dropped  in  the  plane  of  Sq  according 
to  a  distribution  K{Sq)  and  in  such  a  way  that  every  set  Si 
intersects  with  Sq.  Compute  the  fraction  of  Sq  where  at  least 

^Note  that  for  sensing,  we  do  not  require  that  sensors  are  located  within 
the  field  of  interest.  Instead,  as  shown  in  Figure  1(a),  we  require  that  they 
can  monitor  some  part  of  the  field  of  interest  even  if  they  are  located  outside 
of  it. 


y 


Fig.  1.  (a)  A  heterogeneous  sensor  network  with  sensors  covering  the  deployment  region  Ao-  (b)  A  convex  set  A  and  the  corresponding  quantities  that 

define  the  kinematic  density,  (c)  Two  convex  sets  Ao ,  Ai  intersecting,  and  the  common  area  • 


k  out  of  the  N  sets  Si  intersect. 

In  the  mapping  of  the  stochastic  coverage  problem  to  the 
set  intersection  problem,  the  fixed  closed  set  Sq  corresponds  to 
the  field  of  intersect  Aq.  The  N  closed  sets  dropped  according 
to  the  distribution  K{Sq)  correspond  to  the  sensing  areas  of 
the  N  sensors  deployed  according  to  the  distribution  K{Aq). 
By  computing  the  fraction  of  the  set  5o,  where  at  least  k  out 
of  N  sets  Si  intersect,  we  equivalently  compute  the  fraction 
of  the  field  of  interest  that  is  k-covered^. 

The  set  intersection  problem  has  been  a  topic  of  research  of 
Integral  Geometry  and  Geometric  Probability  [7],  [14],  [17]- 
[19].  Before  we  provide  analytical  coverage  expressions  based 
on  our  formulation,  we  present  relevant  background. 

B.  Background  on  Integral  Geometry 

In  this  section,  we  present  relevant  background  on  Integral 
Geometry  that  we  use  in  Section  IV  for  deriving  analytical 
coverage  expressions  based  on  our  formulation.  Interested 
reader  is  referred  to  [7],  [14],  [17]-[19],  as  reference  to 
Integral  Geometry. 

We  first  define  the  notion  of  the  kinematic  density  for  the 
group  of  motions  of  a  set  A  in  the  plane,  that  is  used  to  define 
a  measure  that  quantifies  the  possible  positions  of  A^  such  that 
a  specific  event  occurs  [17].  The  kinematic  density  expresses 
the  differential  element  of  motion  of  a  set  in  the  plane  and  is 
defined  as  follows. 

Definition  1:  Kinematic  Density-Let  Ai  denote  the  group 
of  motions  of  a  set  A  in  the  plane.  The  kinematic  density  dA 
for  the  group  of  motions  Ai  in  the  plane  for  the  set  A,  is 
defined  as  the  differential  form: 

dA  =  dx  A  dy  A  d(j),  (1) 

where  A  denotes  the  exterior  product  used  in  exterior  calculus 
[5],  [6],  (x,  y)  denote  the  Cartesian  coordinates,  and  denotes 
the  rotation  angle  of  A  with  respect  to  the  x  axis  of  the 
coordinate  system'*  [17]. 

^Due  to  their  equivalence,  ,4o  and  iSo  as  well  as  the  terms  sensing  area 
and  set  are  used  interchangeably  in  the  rest  of  the  paper. 

'^For  every  set  A,  one  can  randomly  choose  a  reference  point  O,  based  on 
which  all  translations  and  rotation  motions  of  A  are  defined. 


In  figure  1(b),  we  show  a  set  S,  a  randomly  selected 
reference  point  O  G  S,  and  the  axis  of  a  coordinate  system. 
All  rotations  and  translations  for  the  set  S  are  defined  with 
respect  to  the  reference  point  O.  Integrating  the  kinematic 
density  of  a  set  A  over  a  group  of  motions  Ai  in  the  plane, 
yields  a  measure  for  the  set  of  motions  Ai,  which  is  called 
the  kinematic  measure  [17],  defined  below. 

Definition  2:  Kinematic  measure-The  kinematic  measure 
m  of  a  set  of  motions  Ai  in  the  plane  is  defined  by  the  integral 
of  the  kinematic  density  dA  over  Ai  : 


By  measuring  the  motions  of  a  set  in  the  plane,  we  quantify 
the  space  of  all  possible  positions  of  the  set  that  correspond 
to  that  motion.  The  quotient  of  the  measure  of  any  random 
motion  path  Z  over  the  measure  of  all  possible  motions  Ai  in 
the  plane,  yields  the  probability  p{Z)  for  that  random  motion 
path  Z  to  occur: 


PiZ) 


m{Z) 

m(Ai) 


(3) 


The  kinematic  measure  allows  us  to  compute  the  geometric 
probability  for  a  specific  set  configuration  to  occur,  as  depicted 
in  (3).  Equation  (3),  is  used  in  our  formulation  to  derive  the 
fraction  of  the  field  of  interest  covered  by  a  sensor  deployment, 
as  it  is  illustrated  in  the  following  section. 


IV.  Coverage  in  heterogeneous  sensor  networks 


In  this  section,  we  derive  analytical  expressions  for  coverage 
by  analyzing  the  coverage  problem  as  a  set  intersection  prob¬ 
lem.  We  first  illustrate  the  coverage  computation  when  a  single 
sensor  is  randomly  deployed  to  monitor  the  field  of  interest, 
by  studying  the  intersection  of  two  sets  in  the  plane.  We  then 
compute  coverage  when  the  sensor  is  deployed  according  to  a 
distribution  K{Ao).  We  extend  our  expressions  to  the  general 
case  where  N  sensors  are  deployed  at  random.  We  compute 
the  fraction  of  the  field  of  interest  covered  by  exactly  k  sensors 
in  the  case  of  heterogeneous  sensing  areas  and  simplify  the 
formula  when  the  sensors  have  identical  sensing  areas.  Finally, 
we  compute  the  fraction  of  the  field  of  interest  covered  by  at 
least  k  sensors. 


A.  Coverage  Achieved  by  Random  Deployment  of  a  Single 
Sensor 

Let  us  consider  the  simple  case  where  a  single  sensor  si 
is  randomly  deployed  in  such  a  way  that  it  monitors  some 
part  of  the  field  of  interest.  The  achieved  coverage  can  be 
computed  by  considering  the  intersection  of  two  sets  in  the 
plane.  Let  Ao,Ai  denote  two  sets  in  a  plane  with  Ao  being 
fixed,  while  Ai  can  move  freely,  represents  the  field  of 
interest,  while  Ai  represents  the  sensing  area  of  node  si.  The 
average  size  of  the  common  area  Aoi  between  sets  Ao,Ai, 
when  Ai  is  randomly  dropped  in  the  plane,  defines  the  area 
of  ^0)  covered  by  Ai-  Normalizing  Aqi  over  Aq  we  obtain 
the  fraction  fr{Ao)  of  Ao  covered  by  Ai.  In  figure  l)(c),  we 
show  two  sets  A/,Ai  and  the  common  area  between  them. 

To  compute  fr{Ao),  we  randomly  select  a  point  P  of  Ao: 
and  find  the  set  of  all  positions  of  Ai  that  include  P.  Dividing 
the  measure  of  all  the  positions  of  Ai  that  include  P  over  the 
measure  of  all  the  positions  of  Ai  such  that  AoC\Ai  f  0 
yields  the  probability  p{P  G  ^i)  that  the  randomly  selected 
point  P  is  covered  by  Ai  [17],  [18].  Integrating  p{P  G  Ai) 
over  all  P  G  normalizing  over  the  size  of  yields 

fr(Ao)-  The  following  theorem  holds  only  for  convex  sets, 
though  it  can  be  extended  in  the  case  of  non-convex  sets  by 
appropriate  computation  of  the  kinematic  measures  [17],  [18]. 


Theorem  1:  Let  be  a  fixed  convex  set  of  area  Fq  and 
perimeter  Lq,  and  let  Ai  be  a  convex  set  of  area  Fi  and 
perimeter  Li,  randomly  dropped  in  the  plane  in  such  a  way 
that  it  intersects  with  Aq.  The  probability  that  a  randomly 
selected  point  P  G  Aq  is  covered  by  Ai  is  given  by: 


p{P  G  Ai) 


2ttFi 

27r(Lo  +  Fi)  +  LqLi 


(4) 


-4o  n  A  7^  0  is: 


dAi 


IAoH-AiAO 


dx  A  dy  A  df 


JAon-Ai^D 

27r(Fo  +  Fi)  +  LqLi.  (6) 


Due  to  the  length  and  complexity,  the  proof  of  (6)  is  omitted. 
Interested  reader  is  referred  to  [17],  [18]  for  details. 

By  combining  (5)  and  (6)  we  can  compute  the  probability 
p{P  G  Ai)  as: 


p{P  G  Ai) 


m(^i  :  F  G  A  n  -4i) 
m{Ai  ■.  AQ[]Aif%) 

27rFi 

27r(Fo  +  Fi)  +  LqLi 


(7) 


Note  that  p{P  €  Ai)  is  only  dependent  on  the  area  and  the 
perimeter  of  the  convex  sets  that  intersect  and  not  on  the  shape 
of  those  sets. 


Lemma  1:  The  fraction  fr{Ao)  of  a  fixed  convex  set  of 
area  Fq  and  perimeter  Lq  that  is  covered  by  a  convex  set  Ai 
of  area  Fi  and  perimeter  Li,  when  Ai  is  randomly  dropped 
in  the  plane  in  such  a  way  that  it  intersects  with  is  given 
by: 


A(A) 


27rFi 

27r(Fo  +  Fi)  +  LqLi 


(8) 


Proof:  Equation  (7)  expresses  the  probability  that  a 
randomly  selected  point  P  G  is  covered  by  Ai .  Integrating 
(7)  over  all  points  P  G  Aq  provides  the  size  Fqi  of  the 
common  area  Aqi  between  A/  and  Ai  : 


Proof:  The  probability  that  P  is  covered  by  Ai  is  equal 
to  the  measure  of  the  set  of  motions  of  Ai  such  that  P  G  Ai 
divided  by  the  measure  of  the  set  of  motions  of  Ai  such  that 
A  n  “4i  f  0.  We  now  compute  the  two  measures. 


m{Ai  :  P  G  Aq  riA) 


w 

f 

dAi 

J  PeAo 

n-4i 

(ii) 

f 

dA\ 

J  PgAi 

r 

p2-K 

= 

/ 

dx  A  dy  /  df 

Jp&Ai 

Jo 

= 

27rFi, 

(5) 

where  in  5(i)  we  integrate  dAi  over  all  motions  of  Ai  such 
that  P  G  -4o  n  .  Since  by  assumption  P  G  Aq  and  is 
fixed,  in  5(ii)  we  integrate  dAi  over  all  motions  of  Ai  such 
that  P  G  Ax-  The  measure  of  all  motions  of  Ai  such  that 


II 

[  p{P  G  Ai)dP 

J  P^Aq 

= 

p{PgAi)  [  dP 

JpeAo 

= 

p{P  G  Ai)Fo 

2'kFqFi 

2'k{Fq  +  Fi)  +  LqLi 

Normalizing  Fqi  by  Fq 

yields: 

A(A)  = 

Pqi 

27rFoFi  1 

27r(Fo  +  Fi)  +  LqLi  Fq 

2'kFi 

27r(Fo  +  Fi)  +  LqLi 

= 

p{PgAi). 

(9) 


(10) 


B.  Coverage  Achieved  by  Deployment  of  a  Single  Sensor 
According  to  an  Arbitrary  Distribution 

In  the  case  where  the  Ai  is  not  randomly  deployed  in 
the  plane,  but  it  follows  an  arbitrary  distribution  K{Aq),  the 
measures  in  (5),  (6)  are  calculated  as  weighted  functions  of 
the  probability  density  function  k{x,y,(l))  of  Ai- 


m{Ai-.A^{^Aif%) 
m{Ai  :  P  G  Ao  riA) 


kdxAdyAdf,  (11) 
kdx  A  dy  A  d(j),(12) 


iPeAi 


where  Z  =  AnC\Ai  f  0.  Depending  on  the  distribution 
K{Af),  the  measures  in  (11),  (12)  may  have  a  closed  form. 
When  ^oo  is  deployed  according  to  the  distribution  K{Ao), 
we  can  calculate  the  probability  p{P  G  ^i),  by  substituting 
the  measures  in  (11),  (12)  into  (7).  The  p{P  G  Ai),  is  the 
basic  building  block  for  deriving  expressions  for  coverage  in 
the  general  case  where  N  sensors  are  deployed,  as  we  show 
in  the  following  section. 


C.  Coverage  in  the  Case  of  Multiple  Sensors 

In  this  section,  we  compute  the  probability  p{S  =  k)  that 
a  randomly  selected  point  P  G  Ao  is  covered  by  k  sensors 
when  N  sensors  are  randomly  deployed.  Using  p{S  =  fc),  we 
compute  the  probability  that  P  is  covered  by  at  least  k  sensors, 
as  well  as  the  fraction  of  covered  by  at  least  k  sensors. 


Theorem  2:  Let  Ao  be  the  field  of  interest  of  size  Fq  and 
perimeter  Lq,  and  let  N  sensors  with  sensing  area  Ai  of  size 
Fi  and  perimeter  Li  be  deployed  over  ^o-  The  probability 
p{S  =  k)  that  a  randomly  chosen  point  P  of  is  covered 
by  exactly  k  sensors  when  fc  >  1  is  given  by: 


p{S  =  k) 


n^Li(27r(To  +  Fr)  +  LoLr) 


(13) 


where  J{i,j)  =  (27rLo  +  LgLci  ^),  T  is  a  matrix  in  which 
each  row  j  is  a  k-permutation  of  [1 . . .  iV],  and  G  is  a  matrix 
in  which  each  row  j  contains  the  elements  of  [1 . . .  A^],  that 
do  not  appear  in  the  row  of  T. 


Proof:  In  order  to  prove  Theorem  2,  we  map  the  problem 
of  coverage  to  the  set  intersection  problem,  as  illustrated  in  our 
problem  formulation  in  Section  III-A.  When  a  single  sensor  Si 
is  deployed,  the  probability  that  it  covers  a  randomly  selected 
point  P  e  ^0  is  given  by  Theorem  1.  Hence,  the  probability 
p(P  f.  Ai)  can  be  computed  as: 

p{PiAi)  =  i-p(pgA) 

_  ^ _ ^TlPj _ 

27r(To  +  Ff)  +  L^Li 
_  27rfo  +  LpAj 

27r(Fo  +  Ff)  + 

Given  that  fact  that  the  N  sensors  are  independently  deployed 
in  the  plane  so  that  they  cover  some  part  of  ?  the  probability 
p{S  =  k)  that  a  randomly  selected  point  P  G  Aq  is  covered  by 


exactly  k  sensors  is  equal  to  the  probability  that  P  is  covered 
by  exactly  k  specific  sets.  Let  T  denote  a  fcx(^)  matrix  where 
each  row  j  is  a  k-permutation  of  the  vector  [1 ...  At],  and  let 
G  denote  a  (At  — fc-|-l)x('^)  matrix  where  each  row  j  contains 
the  elements  of  [1 ...  At],  that  do  not  appear  in  the  row 
of  T.  Consider  for  example,  r(l)  =  [1 . . .  A:]  and  G(l)  = 
[A:  -f  1 ...  At].  The  probability  p(r(l))  that  P  is  covered  by 
exactly  the  sets  with  indexes  in  the  first  row  of  T  is  given  by: 

p(T(l))  i  p{PGAi,...,PiAk+i,...,PiAN) 
p{PGAi),...,{PGAk) 
p{P  Ak+i),  ■  ■  ■  ,p{P  An) 

(iii)  27r  Fi 

2tt{Fq  +  Fi)  +  LqLi 

_ 27rJtfc _ 

2t:{Fq  +  Fk)  +  Lo-Tfc 
2ttFq  +  LoLk+i 
2tt{Fo  +  Fk+i)  +  LoLfc+i 
27rFo  +  LqLn 
2tt{Fo  +  Fn)  +  LqLn 
ati(2trf,  )  nlfc+i(2^T-o  +  LqL,) 
n!Li  (2tr(i^o  +  Fr)  +  LoLr) 

UU  (2ttTT,.,)  i‘27rFo  +  LqLg,  J 

nf=i  (27r(Fo  +  Fr)  +  LoLr) 

(15) 


In  (i),  we  show  which  k  sets  include  point  P.  Due  to  the 
independence  in  the  set  deployment,  in  (ii),  the  intersection  of 
the  events  in  (i)  becomes  a  product  of  the  individual  events.  In 
(iii),  we  substitute  the  individual  probabilities  from  (7),  (14). 
In  the  general  case,  the  probability  that  the  sets  with  indexes 
of  the  row  of  T  cover  point  P  is  given  by: 


p{T{i)) 


nf=i  (MFo  +  Fr)  +  LoLr)  ■ 


(16) 


Since  we  are  not  interested  in  a  specific  set  permutation  to 
cover  point  P,  the  probability  that  p{S  =  k)  is  a  summation 
of  p{T{i))  for  all  possible  A:-permutations.  Summing  p{T{i)) 
over  all  i  yields  (13): 


(”) 

p{S  =  k)  =  ^p(T(z)) 

h  \  nil  (MFo  +  Fr)  +  LoLr)  J 


According  to  Lemma  1,  (13)  also  expresses  the  fraction  of 
that  is  covered  by  exactly  k  sensors.  Equation  (13)  is  valid  for 
A:  >  1.  The  fraction  of  the  Ao  that  is  not  covered  by  any 
sensor,  is  given  by  the  following  corollary. 


Non-covered  fraction  of 


Fig.  2.  Fraction  fr{Ao)  of  Ao,  that  remains  non-covered  as  a  function  of  the  number  of  sensors  N  that  are  deployed  to  monitor  the  field  of  interest. 


Corollary  1:  The  fraction  of  that  is  not  covered  by  any 
sensor  when  N  sensors  are  randomly  deployed  is  given  by, 


p{S 


Q-,  _  TT  27rTo  -p  L^Li 

’  ~  27r(Fo  +  F,)  +  Loi^/ 


(17) 


Proof:  Given  that  fact  that  the  N  sensors  are  indepen¬ 
dently  deployed  in  the  plane  so  that  they  cover  some  part  of 
„4o,  the  probability  p{S  =  0)  that  none  of  the  Ai,  i  =  1 . .  .N 
covers  point  P  is: 


p{S  =  0) 


(i) 


(ii) 


p{P  (ji  -4i, . . . ,  P  ^  An) 

N 

-pr  /  ^TtFq  +  LqLi  \ 
y  \2Tr{Fo  +  F,)  +  LolJ  ' 


(18) 


Equality  in  (i)  holds  due  to  the  independence  in  the  deploy¬ 
ment  of  the  sensors  Si.  In  (ii),  we  substitute  p{P  ^  Ai)  from 
(14).  ■ 


Theorem  3:  Let  Aq  be  the  field  of  interest  of  size  Fq  and 
perimeter  Lq,  and  let  N  sensors  with  sensing  area  Ai  of  size 
Fi  and  perimeter  Li  be  deployed  over  Aq.  The  probability 
that  a  randomly  selected  point  of  Aq  is  covered  by  at  least  k 
sensors  is  given  by: 

fc-i 

p{S>k)  =  l-'^p{S  =  h)  (20) 

h^l 

where 

p{S  =  h)  =  . —  \  — 

nr=l(2^(^0  +  Ff)  +  LqLA) 

Proof:  Theorem  3,  holds  by  observing: 

fc-i 

p{S  >k)  =  1  —  p{S  <  fc)  =  1  —  '^^p{S  =  h),  (21) 

h=l 

and  substituting  (13)  into  (21).  ■ 


In  the  case  where  the  sensors  have  identical  sensing  area,  that 
is,  Fi  —  F  and  Li  =  L  then  the  following  corollary  holds. 


Corollary  2:  Let  Fi  =  F  and  Li  =  L.  The  probability 
that  a  randomly  selected  point  of  Aq  is  covered  by  exactly  k 
sensors  is  given  by 


p{S  =  k) 


(^)(2^P)fe(2^Po  +  Po£)^-fe 
(27r(Po  +  P)  +  LqL)^ 


(19) 


V.  Validation  of  the  Theoretical  Results 

In  this  section,  we  validate  our  theoretical  results  derived  in 
Section  IV  via  simulation.  We  perform  experiments  for  both 
homogeneous  and  heterogeneous  sensor  networks  and  show 
that  the  theoretical  formulas  match  the  simulations.  We  also 
provide  an  example  for  analytically  computing  the  number  of 
sensors  that  need  to  be  deployed  in  order  to  achieve  the  desired 
degree  of  coverage. 


Proof:  Corollary  2  holds  by  substituting  Fi  =  F  and 
Li  =  L,  into  (13).  ■ 

Once  we  have  computed  the  probability  for  a  randomly 
selected  point  P  of  Aq  to  be  covered  by  exactly  k  sensors, 
we  can  also  compute  the  probability  that  a  randomly  selected 
point  P  is  covered  by  at  least  k  sensors. 


A.  Coverage  in  Homogeneous  Sensor  Networks 

In  our  hrst  experiment,  we  randomly  deployed  a  variable 
number  of  sensors  with  identical  sensing  area  in  a  disk  of 
radius  R  =  100m.  All  sensors  had  a  circular  sensing  area  of 
radius  r  =  10m.  We  repeated  the  experiments  100  times  and 
averaged  the  results.  We  first  compute  the  fraction  /r(^o)  of 
Aq,  that  remains  non-covered  as  a  function  of  the  number  of 
sensors  N  that  are  deployed  to  monitor  the  held  of  interest. 


Probability  density  function  p(S=k)  for  N=200 


(a) 


Probability  density  function  p(S=k)  for  N=600 


(C) 


Probability  density  function  p(S=k)  for  N=1 000 


(e) 


Fraction  of  covered  by  at  least  k  sensors  for  N=200 


(b) 


Fraction  of  covered  by  at  least  k  sensors  for  N=600 


(d) 


Fraction  of  A^  covered  by  at  least  k  sensors  for  N=1000 


(f) 


Fig.  3.  (a)  The  pdf  of  the  fraction  fr(Ao)  covered  by  exactly  k  sensors  when  N  =  300  sensors  with  identical  sensing  area  are  randomly  deployed,  (b) 
The  fraction  fr{Ao)  covered  by  at  least  k  sensors  when  N  =  300  sensors  with  identical  sensing  area  are  randomly  deployed,  (c)  The  pdf  of  the  fraction 
fr{Ao)  covered  by  exactly  k  sensors  when  N  =  500  sensors  with  identical  sensing  area  are  randomly  deployed,  (d)  The  fraction  fr{Ao)  covered  by  at 
least  k  sensors  when  N  =  500  sensors  with  identical  sensing  area  are  randomly  deployed,  (a)  The  pdf  of  the  fraction  fr{Ao)  covered  by  exactly  k  sensors 
when  N  =  1000  sensors  with  identical  sensing  area  are  randomly  deployed,  (f)  The  fraction  fr{Ao)  covered  by  at  least  k  sensors  when  N  =  1000  sensors 
with  identical  sensing  area  are  randomly  deployed. 


The  theoretical  formula  that  computes  fr{Ao)  is  obtained 
from  Corollary  1  and  is  equal  to: 


fr{An)  =p{S  =  0) 


/  2ttFo  +  LoL  Y 
{2tt{Fo  +  F)  +  LoL  )  ’ 


(22) 


where  Fq  =  =  2'kR,F  =  7rr^,L  =  27rr.  In  figure 

2(a),  we  show  the  fraction  fr{Ao)  of  Ao,  that  remains  non- 
covered  as  a  function  of  the  number  of  sensors  N  that  are 
deployed  to  monitor  the  field  of  interest.  We  observe  that 


the  theoretical  formula  in  (22)  conforms  with  the  simulation 
results.  Since  our  method  does  not  suffer  from  the  border  effect 
problem,  (22)  is  accurate  despite  the  bounded  size  of  the  field 
of  interest. 

In  figure  3(a),  we  show  the  pdf  of  the  fraction  fr{Ao) 
covered  by  exactly  k  sensors  when  N  —  200  sensors  with 
identical  sensing  area  are  randomly  deployed.  The  equivalent 
sensor  density  is  equal  to  p  =  0.0063  sensors/m^.  The 
same  graphs  for  N  =  600,  iV  =  1000  (densities  p  =  0.019 


Fig.  4.  Fraction  fr{Ao)  of  ^o,  that  remains  non-covered  as  a  function  of  the  number  of  sensors  N  that  are  deployed  to  monitor  the  field  of  interest. 


sensors/m^,  p  =  0.032  sensors/m^)  are  provided  in  figures 
3(c)  and  3(e),  respectively.  The  pdf  of  /r(.4o)  is  equal  to 
the  probability  that  a  randomly  selected  point  P  is  covered 
by  exactly  k  sensors.  Our  analytical  derivation  in  Section  IV, 
yields: 

fr{AQ)  =  p{S  =  k) 

(^)(2^f)fe(2^f-o  +  Lo£)^-fe 
{2TT{Fo  +  F)+LoLf 

In  figure  3(b),  we  show  the  fraction  of  covered  by 
at  least  k  sensors  when  N  =  200.  The  same  graphs  for 
N  =  600,  N  =  1000  are  provided  in  figure  3(d)  and 
3(f),  respectively.  For  both  values  of  N  we  observe  that  our 
theoretic  formulas  conform  with  the  simulation  results.  For  all 
graphs  in  figure  2,  3  we  show  the  theoretical  result  according 
to  our  expressions,  and  the  simulation  values. 


B.  Coverage  in  Heterogeneous  Sensor  Networks 

In  our  second  experiment,  we  considered  a  hierarchical 
(heterogeneous)  sensor  network,  where  two  types  of  sensors 
are  deployed.  Type  A  has  a  sensing  area  of  disk  shape  with  a 
sensing  range  ta  =  lOm,  while  type  B  has  a  sensing  area  of 
disk  shape  with  a  sensing  range  of  =  15to.  We  randomly 
deployed  an  equal  number  Na  =  Nb  =  y  of  sensors  of 
each  type  over  a  circular  field  of  interest  of  size  Fq  =  itF? 
where  R  —  lOOm.  In  figure  4,  we  show  the  fraction  fr{Ao) 
of  Ao,  that  remains  non-covered  as  a  function  of  the  number 
of  sensors  N  that  are  deployed  to  monitor  the  field  of  interest. 
The  theoretical  formula  that  compute  that  is  equal  to: 


fr{Ao)  =  p{S  =  0) 

_  TT  +  LqLi 

y  2^(Fo  +  F,)  +  Loi*’ 


(24) 


where  Fq  =  ttR^,  Lq  =  27ri?,  F)  =  Trrf ,  L  =  27rri. 

We  observe  that  the  simulation  results  verify  the  validity 
of  our  theoretical  expression.  In  figure  5(a),  we  show  the  pdf 
of  the  fraction  fr{Ao)  covered  by  exactly  k  sensors  when 
N  =  300  sensors  are  randomly  deployed.  The  equivalent 


sensor  density  is  equal  to  p  =  0.0095  sensors/m^.  The 
same  graph  for  N  =  500,  =  1000  (densities  p  =  0.019 

sensors/m^,  p  =  0.032  sensors/m^)  are  provided  in  figures 
5(c)  and  5(e),  respectively.  The  fr{Ao)  covered  by  exactly  K 
sensors  is  equal  to  the  pdf  p{S  =  k)  of  the  probability  that 
a  randomly  selected  point  P  is  covered  by  exactly  k  sensors. 
Our  analytical  derivation  in  Section  IV,  yields: 

fr{Ao)  =  p{S  =  k) 


For  fc  =  0  : 


27rFn  +  LoLi 


( 

^  n  (^27r(Fo  +  F,)  +  Fo 


while  for  fc  >  1  : 

eEJ  (n,ti(27rFT.,,)nf=i'=^(F^)) 

fr{Ao)  =  - 

nr=,{27r{Fo  +  Fr)  +  LoLr) 

In  figure  5(b),  we  show  the  fraction  of  Aq  covered  by 
at  least  k  sensors  when  N  =  300.  The  same  graphs  for 
N  =  500,1V  =  1000  are  provided  in  hgures  5(d),  and 
5(f),  respectively.  We  again  verify  that  our  theoretical  formula 
agrees  with  the  simulation  results. 

In  the  case  of  heterogeneous  sensor  networks  where  each 
sensor  has  a  different  sensing  area,  the  formula  in  (25) 
has  an  exponentially  increasing  computational  cost,  since  an 
exponentially  increasing  summation  of  terms  must  be  com¬ 
puted  in  order  to  derive  the  exact  coverage  achieved.  Such  a 
computation  may  not  be  feasible  for  large  networks.  In  such 
a  case,  an  approximation  can  be  used  for  our  formulas  by 
employing  the  expressions  derived  for  a  homogeneous  sensor 
network  and  substituting  the  size  F  and  perimeter  L  of  the 
sensing  area  of  the  sensors  with  the  expected  size  F[F]  and 
expected  perimeter  E[L].  The  theoretical  approximation  for 
such  a  case  is: 


fr{Ao)  = 


piS  =  k) 

Q{2nE[F])>^{2nFo- 


LoE[L]) 


N-k 


{2tt{Fo  +  E[F])  +  LoE[L]) 


N 


-.(25) 


Pdf  p(S=k)  for  N=300  and  heterogeneous  network 


(a) 


Pdf  p(S=k)  for  N=500  and  heterogeneous  network 


(C) 


Pdf  p(S=k)  for  N=1000  and  heterogeneous  network 


(e) 


fr{AQ)  covered  by  at  least  k  sensors  for  N=300 


(b) 

fr{AQ)  covered  by  at  least  k  sensors  for  N=500 


(d) 


fr(AQ)  covered  by  at  least  k  sensors  for  N=1000 


(f) 


Fig.  5.  Heterogeneous  sensor  network,  with  the  field  of  interest  being  a  disk  of  radius  R  =  100m.  An  equal  number  of  two  types  of  sensors  are  deployed; 
Type  A  has  a  sensing  area  of  a  disk  shape  with  radius  =  10m,  while  type  B  has  a  sensing  area  of  a  disk  shape  with  rg  =  15m.  (a)  The  pdf  of  the 

fraction  fr{Ao)  covered  by  exactly  k  sensors  when  N  =  300  sensors,  (b)  The  fraction  fr{Ao)  covered  by  at  least  k  sensors  when  N  =  300  sensors,  (c) 
The  pdf  of  the  fraction  fr{Ao)  covered  by  exactly  k  sensors  when  N  =  500  sensors,  (d)  The  fraction  fr{Ao)  covered  by  at  least  k  sensors  when  N  =  500 
sensors,  (e)  The  pdf  of  the  fraction  fr{Ao)  covered  by  exactly  k  sensors  when  N  =  1000  sensors,  (f)  The  fraction  fr{Ao)  covered  by  at  least  k  sensors 
when  N  =  1000  sensors. 


In  figure  6(a)  we  show  the  pdf  obtained  via  simulation  for 
our  heterogeneous  sensor  network  experiment,  for  N  —  500 
sensors,  the  theoretical  values  based  on  the  exact  formula  in 
(25),  and  the  approximation  in  (25).  In  hgure  6(b),  we  show 
the  fraction  of  Aq  covered  by  at  least  k  sensors.  We  observe 
that  for  the  case  of  heterogeneous  sensor  networks  where 
each  sensor  has  a  different  sensing  area,  (25)  provides  a  good 
approximation  of  the  coverage  achieved,  without  incurring  the 
computational  cost  of  (25). 


C.  An  Example  of  Computing  the  Coverage  in  a  Sample 
Network 

In  this  section,  we  provide  an  example  of  applying  our 
results  to  a  sample  sensor  network.  Consider  an  Fol  of  size 
Fq  =  10®m^  and  perimeter  Lq  =  4, 000m  where  sensors  of 
identical  sensing  area  F  —  IOOtt  and  perimeter  L  —  207r 
are  randomly  deployed.  We  want  to  compute  the  number  of 
sensors  needed  in  order  for  a  randomly  selected  point  of  the 
Fol  to  be  covered  by  at  least  one  sensor  with  a  probability 


Pdf  p{S=k)  for  N=500  and  heterogeneous  network 


(a) 


fr(A|^)  covered  by  at  least  k  sensors  for  N=500 


(b) 


Fig.  6.  Heterogeneous  sensor  network,  with  the  field  of  interest  being  a  disk  of  radius  R  =  100m.  An  equal  number  of  two  types  of  sensors  are  deployed; 
Type  A  has  a  sensing  area  of  a  disk  shape  with  radius  =  10m,  while  type  B  has  a  sensing  area  of  a  disk  shape  with  rs  =  15m.  (a)  The  pdf  of  the 

fraction  fr{Ao)  covered  by  exactly  k  sensors  when  N  =  500  sensors,  (b)  The  fraction  fr{Ao)  covered  by  at  least  k  sensors  when  N  =  500  sensors. 


Pc  =  95%.  Or  alternatively,  the  number  of  sensors  N  needed, 
so  that  a  fraction  pc  =  0.95  of  the  held  of  interest  is  covered 
by  at  least  one  sensor.  Corollary  1  yields; 


p{S>l)  =  l-p{S  =  Q) 

—  1  _  TT  (  27rfo  +  LqL  \ 
"  H  1^27r(Fo  +  F)  +  LoL; 

^  /  27rfo  +  Lo£  Y 

\2Tr{Fo  +  F)  +  LoL  ) 


We  want  to  the  probability  of  1-coverage  to  be  at  least  p{S  > 
1)  >  p.  Hence, 


P{S  >  1)  = 
N  > 


1- 


/  ^ttFq  +  LqL 


\^27r(_Fo  H" 
log  (1  -  Pc) 


log  ( 


27rFQ-\-LoL 

27r(Fo+F)+LoT 


N 


>  PC 


Substituting  the  values  for  pc,  Fq,  Lg,  F,  L  yields  N  >  9, 728 
sensors. 


VI.  CONCLUSION 

We  studied  the  problem  of  stochastic  coverage  in  planar 
heterogeneous  sensor  networks.  We  formulated  the  coverage 
problem  as  a  set  intersection  problem  and  used  results  from 
Integral  Geometry  to  obtain  analytical  expressions  for  the 
coverage  achieved  by  the  deployment  of  N  sensors.  Our 
formulation  generalizes  to  a  heterogeneous  sensing  model 
where  each  sensor  has  a  different  sensing  area,  while  it 
does  not  suffer  from  the  border  effects  problem.  Furthermore, 
our  approach  applies  to  sensor  deployment  according  to  any 
distribution.  To  verify  our  results,  we  performed  extensive 
simulation  and  showed  that  the  simulation  conforms  with  our 
theoretic  formulas. 
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